July 19, 2012 0:6 WSPC/INSTRUCTION FILE spatial'ssc 



International Journal of Modern Physics: Conference Series 
© World Scientific Publishing Company 



Of 



A SPATIALLY RESOLVED SSC SHOCK-IN-JET MODEL 



Stephan Richter* 

ITPA, Universitat Wurzburg, Emil- Fischer- Strafie 31 
Wiirzburg, 97074, Germany 

Felix Spaniert 

ITPA, Universitat Wurzburg, Emil- Fischer- Strafie 31 
Wurzburg, 97074, Germany 



Received Day Month Year 
Revised Day Month Year 



In this paper a spatially resolved, fully self-consistent SSC model is presented. The ob- 
servable spectral energy distribution (SED) evolves entirely from a low energetic delta 
distribution of injected electrons by means of the implemented microphysics of the jet. 
These are in particular the properties of the shock and the ambient plasma, which can 
be varied along the jet axis. Hence a large variety of scenarios can be computed, e.g. the 
acceleration of particles via multiple shocks. Two acceleration processes, shock acceler- 
ation and stochastic acceleration, are taken into account. From the resulting electron 
distribution the SED is calculated taking into account synchrotron radiation, inverse 
Compton scattering (full cross section) and synchrotron self absorption. The model can 
explain SEDs where cooling processes are crucial. It can verify high variability results 
from acausal simulations and produce variability not only via injection of particles, but 
due to the presence of multiple shocks. Furthermore a fit of the data, obtained in the 
2010 multi-frequency campaign of Mrk501, is presented. 
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1. Introduction 

Synchrotron Self Compton (SSC) models have been quite successful in explaining 
the emissions of blazars. However, recent observational results have led to the con- 
clusion that the approach of modeling blobs in blazar jets as homogeneous regions 
employing SSC codes is not sufficient. Especially the observation of intra day vari- 
ability introduces strong constraints for both the size of the emission region and 
the Doppler factor that usually lack observational evidence or even contradict those 
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(e.g. Ref. [T]and Ref. [5]). Furthermore the observed time delays between light curves 
of different frequencies, especially soft lagi^, are hardly explainable by such models. 

Two-zone- models like Ref. |4| can produce scenarios that can reproduce observa- 
tional data by merely varying physical parameters. Although this is an enormous 
improvement in comparison to one-zone-models, into which the electron distribution 
enters as a couple of unphysical parameters, there are still problems: 

The stochastic character of the dominant acceleration process (Fermi-I) is not 
taken into account. It is modeled as a convection in momentum space and enters as 
a term in the Fokker-Planck equation. Furthermore the radiation zone is modeled 
homogeneously. Therefore variations, for example introduced via particle injection 
from the acceleration zone, will affect the entire space immediately. This behav- 
ior will lead to shortened timescales. Finally the cooling of the particles happens 
everywhere quasi locally and not depending on distance to the acceleration site. 

Additional to addressing the aforementioned issues spatially resolved models 
can predict the morphology of the emission region. Using VLBI this is in principle 
observable. 

2. Model 

The here presented model connects the dominant acceleration process with the 
geometry of the simulation box, resulting from the spatial discretization. 

2.1. Geometry 

The discretization along the z-coordinate, e.g. along the jet axis, is shown schemat- 
ically in Fig. [1] The resulting slices are indexed with i. For each of them a bulk ve- 



Fig. 1. A schematic illustration of the used geometry. The spatial discretization (index i) is along 
the jet axis. A shock is represented as a jump in velocity of the background plasma. 

locity is calculated depending on the positions and properties of the shocks, namely 
the compression ratio R and shock speed Vs- The downstream velocity behind each 
shock Vp, expressed in the upstream frame is 




V P 



Vs(R-l) 
R-V2 



(1) 
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The Fermi-I acceleration then evolves naturally if particles undergo convection along 
z as well as pitch angle scattering. In order to obtain a modeling for this process 
while keeping the momentum space one dimensional, we integrate over the pitch 
angle in each half-space separately. Hence we obtain two particle distributions - n\ 
for particles moving downstream and n~ for particles moving the opposite direction. 
The pitch angle diffusion is then modeled as scattering from one half-space into the 
other by a rate that enters the model as a parameter. 



2.2. Kinetic Equations 

To derive the kinetic equation all relevant processes have to be included into the 
Fokker-Planck equation, which is then discretized in z. After performing the inte- 
gration as described above one obtains the convection equation in space as well as 
the decoupled kinetic equation 



dt = 



D f . ^pl + tfrf _ 2D 1 + P JC ( 7 )) • n( 7 ) 



S(*,7,t) , (2) 



that describes the time evolution of the electron distributions nf , i.e. for each 
spatial slice and half-space. In Eq. [^acceleration due to Fermi-II is described via the 
momentum diffusion coefficient D. Furthermore losses due to synchrotron emission 
(/3 S ) and inverse compton scattering (Pic) ar e included. S stands for the source 
function. 

To obtain the SED in each slice, the photon density is calculated employing 
dN 4tt , N 

-7T- = -C ■ K U SSA ■ N + 7— ■ { e v,IC + e^sync) ~ 7 ■ (3) 

dt hv t e 



The dominant gain e v sync is due to synchrotron emission and calculated employing 
the Melrose-ApproximatiorM The synchrotron self absorption coefficient k v .ssa is 
calculated with the delta approximation. To obtain the rates for inverse Compton 
scattering e„jc the full Klein-Nishina cross sectiorP is used. The catastrophic loss 
is parameterized by the escape timescale t esc . The full definitions of all used terms 
can be found in e.g. Ref. [71 Finally the total SED is calculated using the model of 
Blandford and KoniglSI, taking into account light travel times and time dilation. 



3. Results 

The Blazar Mrk501, a HBL (high frequency peaked BL-Lac object) with a redshift 
of z = 0.034, was observed in a multi-frequency campaign between March and 
August 2009^1 Our fit of the averaged low state of this source is shown in Fig. [2] 
This fit was realized in a simulation region of 10 15 cm with a magnetic background 
field of 1 G. Particles were injected with a rate of ~ 200 cm' 3 s _1 and a Lorentz 
factor of ji n j = 345 far upstream i.e. at the edge of the simulated region. The 
shock responsible for the acceleration had a velocity v$ = 0.2 c and a compression 
ratio R — 3. The Doppler factor of the shock frame is 6 = 20. Although the size 
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Fig. 2. Fit of the data obtained during a multi-frequency campaign. The comparison to the two- 
zone model fit, provided by Ref. 2] shows much better agreement, especially in the radio regime. 

of the simulation region is not resolvable in observations, analysis of the produced 
morphology show inhomogeneities that affect the overall SED. 

4. Conclusion 

With the here presented model it is possible to confirm results regarding the particle 
acceleration from one- and two-zone-models and therefore trace them back to the 
jets microphysics. The number of parameters is not higher than that of homogenous 
SSC models and could be further reduced by connecting momentum diffusion and 
the scattering parameter, since both processes occur due to the presence of Alfvcn 
waves. In the fit of the Mrk501 data the difference to homogeneous models is most 
prominent in the radio range. This effect is presumably due to the emission of 
electrons far away from the shock that are already cooled. 

Furthermore our model is able to explain very short time variability since the 
lower limit for the timescale, due to the light crossing time of the emission region, 
no longer holds. Fitting of actual lightcurves will be subject to future publications. 
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